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ABSTRACT 

We report the discovery of 7-ray pulsations (> 0.1 GeV) from the young radio and X-ray pul- 
sar PSR J0205+6449 located in the Galactic supernova remnant 3C 58. Data in the 7-ray band 
were acquired by the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope (formerly 
GLAST), while the radio rotational ephemeris used to fold 7-rays was obtained using both the Green 
Bank Telescope and the Lovell telescope at Jodrell Bank. The light curve consists of two peaks sep- 
arated by 0.49 ± 0.01 ± 0.01 cycles which are aligned with the X-ray peaks. The first 7-ray peak 
trails the radio pulse by 0.08 ± 0.01 ± 0.01, while its amplitude decreases with increasing energy as 
for the other 7-ray pulsars. Spectral analysis of the pulsed 7-ray emission suggests a simple power 
law of index —2.1 ± 0.1 ± 0.2 with an exponential cut-off at 3.0"^' \ ± 0.4 GeV. The first uncer- 
tainty is statistical and the second is systematic. The integral 7-ray photon flux above 0.1 GeV is 
(13. 7± 1.4±3.0) x 10~ 8 cm~ 2 s _1 , which implies for a distance of 3.2kpc and assuming a broad fan-like 
beam a luminosity of 8.3 x 10 34 ergss _1 and an efficiency r\ of 0.3%. Finally, we report a 95% upper 
limit on the flux of 1.7 x 10~ 8 cm~ 2 s _1 for off-pulse emission from the object. 
Subject headings: pulsars: general — stars: neutron 
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1. INTRODUCTION 

The radio source 3C 58 was recogni zed early to be a su- 
pernova remnant (SNR G130.7+3.1: lCaswelilll970ft . and 
subsequently cla ssified as a pulsar wind n ebula (PWN, 
also p lerion) by (jWeiler &: Panagialll978ft . iBecker et al.l 
(1982) identified an X-ray point source in the heart of 
3C 58 as a likely pulsar and subsequent stu dies yielded 
a distance of 3.2 kpc (jRoberts et al.l IT993D . The pul- 
sar J0205+6449 was finally discovered in Chandra X-ray 
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Observatory data, with a period of 65.7 ms, while Rossi 
X-ray Timing Explorer archives allowed a measurement 
of th e spindown rate of P — 1.93 x 1CP 13 (|Murrav et al.1 
2002). Both telescopes observed an X-ray profile with 
two narrow peaks separated by 0.5 in phase. This was 
followed by the detection of weak radio pulsations with 
a pulse averaged flux density of ~ 45/zJy at 1.4 GH z 
and a sharp pulse of width 2 ms (|Camilo et al.l [2002) . 
The pulsar has a very high spin-down luminosity of 
2.7 x 10 37 ergs s _1 (the third most energetic of the known 
Galactic pulsars), a surface magnetic field strength of 
3.6 x 10 12 G, and a characteristic age of 5400 years. It 
also e xhibits a high level of timing noise (R ansom et al.l 
|2004|) . and at least two glitches ha ve occurred since its 
dis covery ([Livingstone et al.l 12008ft . Recently, a study 
bv lLivingstone et alj ( 2009ft presented the first measure- 
ment of the phase offset between the radio and X-ray 
pulse, showing that the first X-ray peak lags the radio 
pulse by = 0.10 ±0.01. 

The 3C 58/J0205+6449 system coincides positionally 
with the 8 28-year-old h istoric al su pernova SN 1181 ac- 
cordin g to lStephens"onl (|197lD and IStephenson fc Green! 
(2002). H owever, recent investigations of models for 
the PWN (IChevalierl I2005D and the velocities of both 
the radio expansio n (|Bietenholz 1120061 ) and optical knots 
(|Fesen et al.ll2008h imply an age for 3C 58 of several thou- 
sand years, closer to the characteristic age of the pulsar 
than of SN 1181. 3C 58 is similar to the apparently 
comparably aged Crab nebula (remnant of SN 1054) 
both having a flat-spectrum radio nebula, non-thermal 
extended X-ray emission, and point-like X-ray emission 
due to a central pulsar. But the two objects differ signif- 
icantly both in luminosity and in size. The radio nebula 
of 3C 58, although ~ 2 times larger, is l ess luminous 
than the Crab by an order of magnitude (jlvanov et al.l 
12004ft. while its X -ray luminosity is ~ 2000 times smaller 
( Torii et ai1l2000D . These disparities could be explained 
by a different age. 

2. OBSERVATIONS 

2.1. Radio Timing Observations 

As part of the sp ace and ground ba sed timing program 
supporting Fermi (|Smith et al.|[200l . PSR J0205+6449 
is being observed at the NRAO Green Bank Telescope 
(GBT) and the Lovell telescope at Jodrell Bank. The 
GBT provides higher precision data, but observations 
are made more frequently at Jodrell Bank, and the ro- 
tational ephemeris used here to fold gamma-ray pho- 
tons is based on times-of-arrival (TOAs) obtained from 
both telescopes between 2008 June 17 and 2009 March 
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9. There are 17 such TO As from GBT mainly at a cen- 
tral frequency of 2.0 GHz, with average uncertainty of 
0.27ms, each based on a lhr integration. The 51 Jo- 
drell Bank TOAs, with average uncertainty of 0.37ms, 
arc derived from observations at 1.4 GHz lasting typi- 
cally for 2 hr. The timing of PSR J0205+6449 is very 
noisy, and in order to describe its rotation well during 
the 9 month interval we use TEMPO 57 to fit to the 
rotation frequency and its first seven derivatives, with 
an rms of 0.4 ms. The best determination of dispersion 
measure (DM =140.7 ± 0.3pccm -3 ) remains that from 
ICamilo et al.l (|2002l ). which we use to correct 2 GHz ar- 
rival times to infinite frequency, with an uncertainty of 
0.3 ms, for comparison with the gamma-ray profile. Note 
that higher signal-to- noise ratio profiles of the pulsar are 
presented at 0.8 and 1.4 GHz in ICamilo et a l. (2002). 

2.2. Gamma-Ray Observations 

The Large Area Telescope (LAT) on Fermi, launched 
on 2008 June 11, is a pair-conversio n telescope consistin g 
of 16 towers set into a 4 x 4 grid (|Atwood et alj [2009) . 
Each tower consists of both a converter-tracker (direc- 
tion measurement of the incident 7-rays) and a CsI(Tl) 
crystal calorimeter (energy measurement/shower devel- 
opment image). The array is surrounded by segmented 
plastic scintillator (charged-particle background identi- 
fication) and connected to a programmable trigger and 
data acquisition system. The instrument is sensitive to 
photons from 0.02 to 300 GeV over a ~ 2.4 sr field of view 
resulting from the compact height /width ratio of the in- 
strument. The LAT hardware design, event reconstruc- 
tion algorithms, background selections and event qual- 
ity selections determine the instrument performance 58 : 
a large effective area on axis (~ 0.8 m 2 ); superior angu- 
lar resolution (#68 ~ 0.5° at 1 GeV for events in the front 
section of the tracker); and an energy resolution better 
than 10% between 0.1 and 10 GeV on axis. The software 
timing chain deriving from a GPS clock on the satellite 
and the phase-folding software h ave been shown t o be 
accurate to better than a few fis ([Smith, et al.ll2008h . 

In this letter, the data collected for the timing analysis 
were obtained by the Fermi LAT in two different ob- 
serving modes, from 2008 July 15 to July 29 during the 
six-week calibration phase, and from 2008 August 3 to 
2009 March 9 when LAT was operating in scanning mode 
under nominal configuration. For the spectral analysis, 
only data acquired from 2008 August 3 were used to avoid 
mixing the different configurations. We used the "Dif- 
fuse" class events having the highest probability of being 
photons. In addition, we excluded the events with zenith 
angles greater than 105° due to the Earth's bright 7-ray 
albedo, and time periods where the Earth's limb came 
within 28° of the source. 

3. ANALYSIS AND RESULTS 

The events were analyzed using the standard soft- 
ware package Science Tools (ST) for the Fermi LAT 
data analysis 59 and the pulsar timing software TEMP02 
(|Hobbs et alj |2006). The timing parameters used in this 

57 http:/ /www. atnf.csiro.au/research/pulsar/tempo/ 



work will be made available on the servers of the Fermi 
Science Support Center 60 . 

3.1. Pulsed Light Curve 

The pulsar is located in the Galactic plane where 
the diffuse gamma radiation is intense, and 5.3° from 
the bright LAT source 0FGL J024 0.3+6113 coinciden t 
with the X-ray binary LSI+61°303 (jAbdo et al.ll2009al ). 
For the timing analysis, a set of photons with energies 
over 0.1 GeV was selected within an energy-dependent 
cone of radius #68 ^ 0.8 x E^y 5 degrees, but with a 
maximum radius of 1.5° with respect to the X-ray pul- 
sar position (I = 130.719°, b = 3.085°). This choice 
takes into account the instrument performance and maxi- 
mizes the signal-to-noise ratio over a broad energy range. 
This truncates the point spread function at low ener- 
gies and decreases the number of background events 
([Atwood et al.ll2009h . A total of 2922 7-rays remain after 
these cuts. The arrival times of events were corrected to 
the Solar System B arycenter using the JPL DE405 Solar 
System ephemeris (Standish ll998l ). and the events have 
been folded using the radio ephemeris from GBT and the 
Lovell telescope at Jodrell Bank. 

Figure [T] (top panel) shows the 50 bin 7-ray phase 
histogram on which we fit each peak. The first peak 
(PI) is offset from the radio pulse (bottom panel) by 
0.08 ± 0.01 ± 0.01 according to a Lorentzian fit with a 
full width at half maximum (FWHM) of 0.15 ±0.03. The 
first phase uncertainty arises from the 7-ray fit, and the 
second is from the DM uncertainty in extrapolating the 
radio time-of-arrival to infinite frequency. The second 
peak (P2) is asymmetric, fit by a two-sided Lorentzian to 
take into account the different widths for the leading and 
trailing edges. The fit places the peak at 0.57±0.01±0.01 
with an FWHM of 0.13 ± 0.04. The two peaks are sep- 
arated by 0.49 ± 0.01 ± 0.01 in phase. We defined the 
'off pulse region' as the pulse minimum between 0.65 
and 1.0 in phase. The dashed line in Figure [1] repre- 
sents the background counts measured from a 2 - 3° ring 
surrounding the pulsar during the pulse minimum. Ev- 
idence for an excess between PI and P2 appears with a 
significance (signal / ^/background) of 5.0 cr between 0.14 
and 0.46 in phase. However, the data cannot constrain 
the pulse phase dependence of this excess flux. 

To examine the energy-dependent trend of the 7-ray 
pulse profile, 50 bin phase histograms are plotted in three 
energy intervals, 0.1 - 0.3 GeV, 0.3 - 1 GeV, and >1 GeV 
in Figure [T] (middle panels), which show an evolution in 
the shape. Between 0.1 and 1 GeV, we see distinctly two 
peaks and a possible bridge region in excess of the un- 
pulsed emission, while above 1 GeV only P2 is significant. 
We also observe that the ratio P1/P2 (sum of counts) de- 
creases with increasing energy, with a ratio of 0.63 ±0.06 
between 0.1 - 0.3 GeV, 0.55 ±0.05 b etween 0.3 -1.0 an d 
0.24 ± 0.06 above 1 GeV. The Vela (lAbdo et alJl2009bh . 
Crab, Geminga and B19 51±32 pulsars show similar be- 
havior ([Thompsonl 1200 ll ) . This general trend suggests a 
spectral energy dependence of the 7-ray light curve. Fi- 
nally, we note that the highest energy photon is in P2 
with an energy of 8.6 GeV. 

Figure [1] (bottom panels) also shows the 2 
60 keV X-ray phase histogram measured by RXTE 



58 http: / / www— glast.slac.stanford.edu / software /IS / glast_lat_performance.htm 

59 http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/overview. html 60 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
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Fig. 1. — Top panel: Phase-aligned histogram of PSR J0205+6449 above 0.1 GeV and within an energy-dependent circle. Two rotations 
are plotted with 50 bins per period. The dashed line shows the background level, as estimated from a ring surrounding the pulsar during 
the off-pulse phase (46 counts/bin). Three following panels: Energy dependent phase histograms for PSR J0205+6449 in the three 
indicated energy ranges , each displayed with 50 bins per pulse period. Second panel from bottom: Count rate in the energy band 2—60 
keV from RXTE data (Livingstone et al. 2009). Bottom panel: Radio pulse profile based on 3.8 hours of GBT observations at a center 
frequency of 2 GHz with 64 phase bins. 
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(|Livingstone et aLll2009f ) as well as the 2 GHz radio pro- 
file from GBT (the sing le peak defines phase = 0). 
iLivingstone et al.l (|2009f ) find that the radio pulse leads 
the X-ray pulse by = 0.10 ± 0.01 and note a separa- 
tion between the two narrow X-ray peaks of A(f> = 0.5. 
The good alignment between the X-ray and 7-ray profiles 
suggests a common origin between the two components. 
This feature is also observed in the LAT data for the Vela 
pulsar where the mor e intense X-ray puls e is aligned with 
the first 7-ray pulse (jAbdo et al.ll2009bf) . 

3.2. Spectra and phase-averaged flux 

To study the phase-averaged spectrum o f J0205+6449, 
a ma ximum likelihood spectral analysis 61 (jMattox et al.l 
1996), implemented in the ST as the gtlike task, was 
performed using a 20° 7-ray map centered on the pulsar 
position between 0.1 and 200 GeV. The systematic errors 
on the effective area are currently estimated as < 5% near 
1 GeV, 10% below 0.1 GeV and 20% above 10 GeV. The 
diffuse emission from the Galactic plane was modeled 
using maps based on the GALPROP model (giLiem.vOl) 
(|Strong et alj l2004aH b(). The extragalactic radiation as 
well as the instrumental backgrounds were modeled by 
an isotropic component with a power-law spectral shape. 
Nearby LAT sources were included in the analysis. For 
PSR J0205+6449, we modeled the shape of the spectrum 
with a power law with an exponential cut-off. 

We first determined the diffuse background compo- 
nents by selecting off-pulse data. Then, we fitted the on- 
pulse data using the diffuse background spectrum thus 
obtained in order to improve the signal-to-background 
ratio of the pulsar and the nearby sources. The best fit 
result is described by: 



dF 
dE 



= N E- T e- E / E ° cm-V^eV" 1 



(1) 



with E in GeV, the term N Q = (1.4 ± 0.1 ± 0.1) x 
lO^cm^s^GeV -1 , the spectral index T = 2.1 ± 0.1 ± 
0.2, and the cut-off energy E c = 3.0tJ;£ ± 0.4 GeV. 
The errors are the statistical and propagated system- 
atic uncertainties, respectively. We obtain from this fit 
over the range 0.1-200 GeV an integral photon flux of 
(13.7±1.4±3.0) x 10 _8 cm~ 2 s _1 and integral energy flux 
of F E ,obs = (6.7 ± 0.5 ± 1.0) x 10- 11 ergscm- 2 s~ 1 . To 
check the assumption of a cut-off energy in the spectrum, 
we have also fit the same dataset with a simple power- 
law of the form dF j dE = N (E/lGcV)- r . The spectral 
model using an exponential cut-off is better constrained 
with a difference between the log likelihoods of ~ 4.5 tr, 
disfavoring the power-law hypothesis. The source, hav- 
ing a statistical significance of 9.5 a (< 10 a) from 0.2 to 
100 GeV for the first three months of the sky survey in 
the science phase of the mission, do es not appear in th e 
Fermi LAT bright 7-ray source list ()Abdo et all 12009 a). 

As a first search for unpulsed emission from the nebula 
3C 58, we fitted a point-source to the off-pulse data at 
the radio pulsar position in the energy band 0.2 - 200 
GeV. No signal was observed from the PWN. Lastly, after 
scaling to the full pulse phase, we derived a 95% CL 
upper limit on the flux of 1.7 x 10~ 8 cm~ 2 s _1 . 

4. DISCUSSION 
61 http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/ 



4.1. Light curve 

Multi-wavelength light curves are important to locate 
the pulse emission in the open field line region and 
hence to understand the mechanisms of particle accel- 
eration. The 7-ray profile for J0205+6449, covering a 
wide range in phase , is reminiscent of the Vela light curve 
(jAbdo et al.l l2009b). including the 7-ray peak alignment 
with the X-ray pulses. The 7-ray delay of 0.08 cycles 
and the peak separation of ~ 0.5 is becoming a consis- 
tent pattern, as t he first 7-ray peaks for Vela, B1951+32 , 
and J2021+3651 ijHalpern et al.l2008HAbdo et al.ll2009cT ) 
lag the radio pulses by 0.13, 0.16, and 0.17 respectively, 
and the separation of the 7-ray peaks is 0.4 - 0.5. This 
fits the predictions of the outer magnetospheric models 
quite well, w hether they be the tradition al outer gap 
model (OG) (g omani fc Yadigaroglul Il995f) or the two 
pole caustic gap model (TPC) (jDvks fc Rudakil2003[ ). 

4.2. Beaming and Luminosity 

To know the 7-ray efficiency, we need to determine 
the total luminosity L 7 = 47r/n(a,C_E) FE,obsD 2 , where 
D is the distance to the pulsar, Fs^obs is the observed 
phase-averaged energy flux for the Earth line of sight 
(at angle C,e to the rotation axis), and /n(a,C-E) is a 
correction factor th at takes into accou nt the beaming 
geometry, given by (|Watters et al.ll200ct ): 



in (a, Cf 



J ' F 1 (a;C,(j))sin(C)dC,d(j) 
2/F 7 (a;C S , 



(2) 



where F 1 {a\C, 1 4 l ) is the radiated flux as a function of 
the viewing angle (Q) and the pulsar phase (0), while 
fn is the ratio between the overall 7-ray emission over 
the full sky and the expected phase-averaged flux for the 
light curve seen from the Earth. In the case of a polar 
cap model ((Ha rding 2005) where charged particles are 
accelerated in charge-depleted zones near the poles of 
the pulsar, /n < 1, while for the outer magnetosphere 
models the emission is radiated with /q > 1. 

An important uncertainty when evaluating the 7-ray 
efficiency arises from the determination of the distance. 
Observations o f neu tral hydrogen (H I) absorption by 
iRoberts et al.1 (|1993l ) yielded an estimate for the radial 
velocity of 3C 58. They convert this to a kinematical dis- 
tan ce of 3.2 kp c assum ing a flat Galactic rotation curve as 
per lFich et alj (|1989f ). who quotes distance uncertainties 
of order 25%, making this r esult consistent with the 2.6 
kpc previously reported by iGreen fc Gulll (fl982h . The 
distance derived fro m the DM is 4.5 kpc, according to 
the NE 2001 model (jCordes fc Lazioll2002D . with an un- 
certainty that can exceed 50%, depending on the viewing 
direction, consistent with the kinematical distance. 

We estimate i 7 = 8.3 x 10 34 (L>/3.2kpc) 2 /n ergss" 1 , 
and deduce, for a neutron star moment of iner- 
tia of 10 45 gcm 2 , a 7-ray efficiency 77 = L 1 /E — 
0.003 /n (-D/3.2 kpc) 2 . Assuming the outer magneto- 
spheric models for PSR J0205+6449, we can deduce 
the parameters a and £ from the 7-ray p eak separation 
and u sing the 7-ray light curve "Atlas" of IWatters et al.1 
(2008). For the OG model, we estimate /n « 0.9- 
1.0 with a - 60°-90° and < - 80°-85°, while for the 
TPC model, f Q w 0.95-1.25 with both a - 50°-90°, 
C - 85°-90° and a - 85°-90°, < - 45°-90°. Examina- 
tion of the geometry of the PWN 3C 58 using Chandra 
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(|Ng fe Romanil l200i [200l yielded a viewing angle of 
C = 91.6 ± 0.2 ± 2.5 ° (inner torus) based on the tilt 
angle of the torus to the plane of the sky. This value is 
in agreement with the results from the OG model and 
more consistent with the first estimated range for the 
TPC model. Adopting an average value for /n = 1 ac- 
cording to the OG model and a distance of 3.2 kpc, we 
evaluate a 7-ray efficiency rj — 0.3% for converting its 
rotational energy loss into 7-rays. 

Studies of rotation powered pulsars suggest that the 
efficiency r\ for converting its rotational energy loss into 
7-rays increases as the open field line voltage V ^2 
4 x io 20 -P~ 3 / 2 P 1 / 2 volts decreases and is proportional 
to E- 1 / 2 (lArond 1 19961). V is al so proportional to the 
open field current dHardind 1 19811) and to the character- 
istic age of the pulsar (|Buccheri et al.l Il978h . We note 
that with an open field line voltage of 1.0 x 10 16 V, 
that PSR J0205+6449 closely follows the approximate 
77 cc E~ x l 2 oc 1/V relation s, confirming the trend. 

iThompson et al.l (|1999l ) point out that for the 7-ray 
EGRET pulsars the broadband energy spectra for the 
young pulsars like Crab and B 1509— 58 peaks in the 
X-ray band, while the older 7-ray pulsars have their 
maximum luminosity in the high-energy regime. Con- 
sidering a X-ray lu minosity between 0.5 and 8 keV of 
1.51 x lO^crgss" 1 (|Kargaltsev fc Pavlovi [20081) and a 
gamma-ray luminosity of 8.3 x 10 34 ergss -1 with /q = 1, 
PSR J0205+6449 seems to have an Lx/Lj ratio closer 
to the middle-aged Vela pulsar than the young Crab and 
B1509— 58 pulsars, and hence suggests that the associa- 
tion of 3C 58 and historical supernova SN 1181 is prob- 
ably incorrect. 

5. SUMMARY 

Using a rotational ephemeris derived from radio ob- 
servations with the GBT and at Jodrell Bank, and 
7-ray data from the Fermi LAT, we have discovered 
in the Galactic SNR 3C 58 7-ray pulsations from 
PSR J0205+6449, the third most energetic of the known 
Galactic pulsars. This source has no EGRET counter- 
part and presents a new opportunity to study the pul- 
sar/SNR association at high energies, demonstrating the 
good performance of the instrument. 



1. The 7-ray profile for J0205+6449 is similar to the 
majority of known 7-ray pulsar light curves, in particular 
to the Vela pulsar. It consists of two peaks separated by 
0.49 ± 0.01 ± 0.01 in phase. The radio pulse leads the 7- 
ray pulse by (j> — 0.08 ± 0.01 ± 0.01 and the X-ray peaks 
are aligned with the 7-ray peaks. 

2. We report a 95% CL upper limit on the flux of 1.7 
x 10 _8 cm _2 s _1 for a possible unpulsed emission from 
the nebula 3C 58. 

3. The 7-ray energy spectrum above 0.1 GeV can be 
described by a simple power law with photon index of 
2.1 ± 0. 1 ± 0.2 and an energy cut-off of 3.01^ ± 0.4 GeV. 

4. Adopting a distance of 3.2 kpc and assuming a broad 
fan-like beam, we report an efficiency of ~0.3% for the 
conversion of spin-down energy into 7-ray emission. 
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